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ABSTRACT 

We discuss the nature of the possible high-£ excess in the Cosmic Microwave Background (CMB) 
anisotropy power spectrum observed by the Cosmic Background Imager (CBI). We probe the angular 
structure of the excess in the CBI deep fields and investigate whether it could be due to the scattering 
of CMB photons by hot electrons within clusters, the Sunyaev-Zeldovich (SZ) effect. We estimate the 
density fluctuation parameters for amplitude, as, and shape, T, from CMB primary anisotropy data 
and other cosmological data. We use the results of two separate hydrodynamical codes for ACDM 
cosmologies, consistent with the allowed as and T values, to quantify the expected contribution from 
the SZ effect to the bandpowers of the CBI experiment and pass simulated SZ effect maps through 
our CBI analysis pipeline. The result is very sensitive to the value of <r 8 , and is roughly consistent 
with the observed power if a$ « 1. We conclude that the CBI anomaly could be a result of the SZ 
effect for the class of ACDM concordance models if as is in the upper range of values allowed by 
current CMB and Large Scale Structure (LSS) data. 

Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

The Cosmic Background Imager (CBI) provides some 
of the highest resolution observations of the cosmic mi- 
crowave background that have been made thus far. The 
observations cover the multipole range 400 < £ < 4000 
which corresponds to collapsed masses ranging from a 
factor ten larger than the local group to the largest super- 
clusters. These observations show, for the first time, the 
fluctuations on scales which gave rise to gala xy clusters 
and the damping of the power on small scales {Silk 1968; 
iPeebles fc Yulll970t iBond fc Efstathioulll987l) . Together 
with the res ults of other recen t , high precision, CMB 
experiments (iMiller et al.1 Il999t iNette rfield et al. 2002; 



|Lee et al.| |2001t |Halverson et al.| [2QQ1D the observations 
fit the scenario of adiabatic fluctuations generated by a 
period of inflationary expansion. The CBI observations 
also provide a unique insight into angular scales where 
secondary anisotropy effects are thought to become an 
important contribution to the power spectrum. 

In a series of papers we have presented the results from 
the first year of CMB observations carried out by the 
CBI between January and December 2000. A prelimi- 
nary analysis of the data was presented in iPadin et al.l 
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(200ll hereafter Paper I). In iMason etnl\ (|2003l here- 
after Paper II) the observations of three differenced 45' 
(FWHM) deep fields were used to measure the power 
spectrum to multipol es £ < 3500 in wide bands A£ ~ 600. 
iPearson et al.l (|2003l hereafter Paper III) discussed the 
analysis of three sets of differenced 145' x 165' mosaiced 
fields. Mosaicing gives a telescope a larger effective pri- 
mary beam than that defined by the dish radii. This in- 
creases the resolution in the iw-plane due to the smaller 
width of the convolving function. It also reduces the ef- 
fect of cosmic variance on the errors by increasing the 
sampled area. The mosaic fields provide high signal- 
to-noise ratio measurements of the power spectrum up 
to multipole s £ ~ 1700 with a resolution A£ = 200. 
iMvers et al.l (1200,1 hereafter Paper IV) give a detailed 
description of our correlation analysis and bandpower 
estimation methods. These have enabled us to analyze 
efficiently the large data sets involved in the CBI mea- 
surements. The implications of th e results on c osmo- 
logical parameters are described in iSievers et alJ {2003, 
hereafter Paper V). 

Results from the deep field observations show a fall 
in the power spectrum up to £ ~ 2000 which is consis- 
tent with the damping tail due to the finite width of the 
last scattering surface. Beyond £ ~ 2000 the power ob- 
served is higher than that expected from standard models 
of damped adiabatic perturbations which provide excel- 
lent fits to the data at larger scales. An extensive set of 
tests have been carried out to rule out possible system- 
atic sources of the measured excess (Paper II). Paper III 
also reports results in this ^-range which are consistent 
with those of Paper II. Since the thermal noise levels in 
the CBI mosaic are substantially higher in this regime 
than in the deep field observations, the current discussion 
will focus on the deep field results. 

The deep field spectrum shows the power dropping to 
£ ~ 2000, beyond which power levels are significantly 
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higher than what is expected based on standard models 
for intrinsic CMB anisotropy. Assuming a single bin in 
flat bandpower above £ = 2010 the observed signal is in- 
consistent with zero and best-fit models at the 3.5<r and 
3.1<j levels respectively. The one and two sigma confi- 
dence intervals for this bandpower are 359-624 /iK 2 and 
199-946 /iK 2 , respectively, with a best-fit power level of 
508 /iK 2 . The confidence limits were obtained by explic- 
itly calculating the asymmetric, non- Gaussian likelihood 
of the high-i? bandpower. 

A key consideration in obtaining these results is 
the treatment of the discrete radio source foreground. 
Sources with positions known from low-frequency radio 
surveys were projected out of the data; the brightest 
sources were also measured at 31 GHz with the OVRO 
40-m telescope and subtracted directly from the data. 
The OVRO data allowed a complementary treatment of 
this potentially limiting population, giving us great con- 
fidence in our bright source treatment. We estimated 
the power level due to sources too faint to appear in the 
low frequency radio source catalogs from 30 GHz num- 
ber counts determined from CBI and OVRO data. The 
contribution of these residual sources is 114 /iK 2 — a fac- 
tor of 4.5 below the observed excess — and we estimated 
an uncertainty of ~ 57 /iK 2 in this contribution. The 
reader will find further discussion of both the CBI deep 
field spectrum and the treatment of the radio source fore- 
ground in Paper II. 

Although the significance of the excess power is not 
conclusive it provides tantalizing evidence for the pres- 
ence of secondary contributions to the microwave back- 
ground anisotropics. One of the strongest expected 
secondary signals is the signature of the scattering of 
CMB photons by hot electrons in clusters known as 
the S unyaev-Zeldovich (SZ) effect ( Sunvaev & Zeldovich 
ll97Q[) . The scattering leads to spectral distortions and 
anisotropics in the CMB. At the CBI observing frequen- 
cies the net effect is a decrement in the temperature along 
the line of sight. (However, as described in §BI the lead- 
trail differencing used in the CBI observations results in 
positive and negative SZ signals in equal measure, so this 
cannot be used as a discriminating signature in our data.) 
The SZ effect is expected to dominate over the primary 
anisotropics at scales of a few arcminutes. 

In this Paper we explore the possibility that the excess 
observed in the CBI deep fields may be a signature of the 
SZ effect. In §Elwe describe constraints on the normal- 
ization of the mass fluctuations as and shape parameter 
T from CMB experiments and a number of independent 
surveys. These parameters are critical in determining the 
amplitude of the SZ effect over the scales of interest to 
the CBI results. In §0|we show the predicted power spec- 
tra for the SZ effect for various cosmological models. We 
use numerical methods to estimate the contribution from 
the SZ effect to the CMB power spectrum. In §0]we use 
simulated CBI observations of SZ 'contaminated' CMB 
realizations to investigate the effect of a SZ foreground 
on our bandpower estimation methods. In §0we extend 
an image filtering technique introduced in Paper IV to 
include specific template filters for the SZ effect. Our 
results and conclusions are discussed in § |H1 



2. AMPLITUDE AND SHAPE PARAMETERS FOR LSS 
FROM CMB DATA 

The SZ power spectrum has a strong dependence on 
the amplitude of the density fluctuations. The amplitude 
is usually parameterized by the rms of the (linear) mass 
fluctuations inside 8/i _1 Mpc spheres, ag. We now sum- 
marize the constraints we can obtain on the amplitude 
and shape of the matter power spectrum from CMB data 
and compare them to results from weak lensing surveys 
and cluster abundance data. Our parameter estimation 
pipeline also includes LSS priors which are designed to 
encompass the range in estimations from such experi- 
ments, as described below. 

We use the parameter T to define the sh ape, fol- 
lowing Bardeen, Bond 1 Kaiser, fc Szalavl ([19861 hereafter 
BBKSrand lEfstathiou et alJ (|1992l hereafter EBW). A 
byproduct of linear perturbation calculations used to 
compute Ci in our database is transfer functions for 
density fluctuations, which can be related to LSS ob- 
servables. Various (comoving) wavenumber scales de- 
termined by the transport of the many species of par- 
ticles present in the universe characterize these spec- 
tra. The most important scale for dark matter dom- 
inated universes is /^ q , that of the horizon at red- 
shift Q m /Q er when the density in nonrelativistic mat- 
ter, £7 m a -3 , equals that in relativistic matter, Vt er a~ A . 
This defines T eq : /c^ q = 5 F~ q h~ x Mpc, where T eq = 

^/i^erA 1 - 68 ^)] -172 - For the cases we consider here 
we simply have fixed the relativistic density to corre- 
spond to the photons and three species of very light neu- 
trinos, so T eq — Q m h. 

Certain functional forms for the transfer functions are 
popular: In EBW a form was adopted which fit a spe- 
cific = 0.03 CDM model, but it is more general to 
adopt a fit to the ^ form given in BBKS appro- 
priately corrected for the difference between the temper- 
ature of_theCMB estimated then and known so well 
now (Bond 1996). Although the coefficients of fits to 
detailed models vary with ujb,h,uj m , which in partic- 
ular result in oscillations in the transfer function for 
large uo^/uom^ it turns out that replacing T eq by T = 
r eg exp[-(^ 6 (l + ^- 1 (2/i) 1 / 2 ) - 0.06)] works reasonably 
well, to about 3% over the region most relevant t o LSS 
( Sugi yamalll995t lBondlll996ft . Further, as shown in lBondl 
(1994), replacing r by r e ff = T + (n s - l)/2 takes into 
account the main effect of spectral tilt over this LSS 
wavenumber band. It has also been shown that the T- 
models do fit the APM and 2dFRS data reasonably well 
at this stage. Particularly exciting is the prospect that 
the baryonic oscillations may be seen, but this is not re- 
quired by the data yet. The approximate codification of 
a vast array of models in a single T e ff variable simplifies 
the treatment of LSS in the CMB data enormously. 

The other variable we use to construct the LSS prior is 
the combination crg^^ 56 . Although various fits to clus- 
ter abundances give slightly different exponents than the 
0.56, the factor is consistent with a number of other mea- 
sures and we adopt this form as a representative value. 
These invariably involve the biasing factor b g for the 
galaxies involved. For example, relating the galaxy flow 
field to the galaxy density field inferred from redshift sur- 
veys takes the form [b g as]P g , where (3 g = Vt^ Q /b g is a 
numerical factor whose value depends upon data set and 
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Fig. 1. — One-dimensional projected likelihood functions of as 
calculated for the CMB data with three prior probability restric- 
tions on the cosmological parameters are contrasted with estimates 
from other datasets. The curves shown used ' all-data' for the 
CMB available for Paper V, na mely DMR DASI lHalverso n et alJ 
l200ll). BOOME RANG ff or thelRuhl et all EPOS) cut), MAXIMA 
ILee et al.ll200ll) . VSA (Scott et al. 20q3) and the CBI mosaic 
data for the odd AL=140 binning. The marginalization has been 
performed over seven cosmological variables and all of the rele- 
vant calibration and beam uncertainty variables associated with 
the experiments, seven in this case. Application of the weak-h 
prior is long-dashed blue, of the flat weak-h prior is solid black 
and of the LSS flat weak-h prior is dotted black. Adding TOCO 
(Miller et al. 199^), the Boomerang test flight and 17 other ex- 
periments predating April 1999 (hereafter 'April 99') as well gives 
very similar results. The Bayesian 50% and associated 16% and 
84% error bars are shown as data points in blue for these and 
other priors, in particular those with the stronger HST measure- 
ment of h included and with SN data included. The magenta data 
points with smaller errors are those determine d with ' all-data' to 
March 2003, including WMAP, as described in lBond etaD (2003). 
The original LSS prior was constructed based on the cluster abun- 
dance data dBondl ll996: Bond & Jaffe 1999t ILange et al J 1.2001V 
An SZ estimate of as from Goldstein et al. (2003|) that simulta- 
neously determines amplitudes for a best-fit primordial spectrum 
and an SZ spectru m for our C BI deep field data in conjunction 
with ACBAR iKuo et aHl200l) and BIMA ( Da wson et all 2002) 
data is shown at the top. Estimates of as from cluster abun- 
dance data and from weak lensing data are also shown. These 
results have invariably had more restrictive priors imposed than 
those for the CMB, and so are not always applicable, but the over- 
all level of agreement in the various approaches is encouraging. 
From t op to bottom, the sampl e cluster value s are from Eke et all 
(1996); Carlbert TeTlH Il997h : iFaTTfc Bahcall (199 ^:IPen| dl998cV 
Pierpaoli, Scott & White (2001); Reiprich & Bohringer (2002); 
Seli ak " 00021) : IViana. Nichol fe Liddlel 120011): IBorgani et alJ 
(2001) , then two estimates fro m Pierpaoli et al. (2003), then from 
Schue ckeTeraLl (120031) ; lAllen et al.H2003l) : IVoevodkin fc Vikhlininl 
(2004). From top to bottom, the weak lensing estimates are from 
Hoeks trTeT aL[(|2002l) ; En Waerbe ke"et al l ([2002): Refregie r"eTaT| 
(2002V "lBacon et al J 120021): IJarvis et alJ J2OO3I) : IHamana et all 
(12001) : lBrown et all (l200^V IHevmans et alJT200f) . 

analysis procedure. A combination such as this also en- 
ters into redshift space distortions. A great advantage of 
the weak lensing and cluster abundance results shown in 
the figures is that they are independent of galaxy bias. 
However, for the cluster case, assumptions are needed 
which, as the spread in estimates indicates, lead to un- 
certainties. 

In Figures ^ and [21 we show the effect of the differ- 
ent priors on the distributions of as and ag^^ 56 , respec- 
tively. We also compare these to the weak lensing results 



Fig. 2. — The weak prior we use for LSS in as^^ 56 * s compared 
with estimates from SZ, cluster abundances, weak lensing and from 
the CMB data shown in Fig. [I] appropriately scaled. The blue 
error bars are for 'all-data' to June 2002, and the magenta are for 
'all-data' to March 2003. The prior used in past work was shifted 
up in central value from 0.47 to 0.55, but was otherwise the same. 
We have also considered the LSS(low-crs) case (dashed green), with 
the prior shifted downward to be centered on 0.40 to accommodate 
better the low cluster abundance estimates, with results shown in 
Tables El 

shown in the figures. T he "Red Cluster Survey" results of 
Hoekstra et alJ (I2QQ2D and the "Virm os-Descart Survey" 
results of Ivan Waerbeke et alJ (|2QQ2fl correspond to a ver- 
sion of a weak prior: they marginalize over T e ff in the 
range 0.1 to 0.4, which is largely equivalent for this appli- 
cation to marginalization over n s and oj m . They also 
marginalize over the uncertainty in the mean redshift of 
the lensed galaxies z s from 0.27 to 0.34 in the former case, 
and 0.78 to 1.08 in the latter case. On the other hand, the 
cosm ology wa s kept fixed with Q m = 0.3. Refre eder et alJ 
(2002| and IB aeon et al J ((2002) adopt a more restrictive 
p arameter range. If th e w eak marginalization scheme 
of lHoekstra et al.1 ()2002ft and Ivan Waerbeke et all (2002) 
were used the error bars would increa se in these cases. 
(Recently |van Waerbeke et al.| (|2005|) have improved 
the Virmos-Descart Surve y analysis, which lowers the 
Ivan Waerbe ke et all (|2002f) cr 8 estimate by 12%.) 

Figs. Hand [2 show sample cluster abundance estimates 
o f as which are discus s ed bel ow . 

Ivan Waerbeke et all ([2002D give a weak lensing result 
of r e ff = 0.25 ± 0.13 with marginalization over Q m from 
0.1 to 0.4 and over as. This is not explicitly shown in 
Fig. [HI where we plot the prior probability we adopted 
for r e ff. The APM result is the long standing one used 
to construct the original prior, that T e ff in the 0.15 to 
0.3 ra nge pr ovided a good fit t o the data, e.g. EBW and 
Bond] Jl99l. Recent 2dF RS peacock et al.ll200l and 
SDSS (Szala v et alJ 12001) results shown give compati- 
ble results. iDodelson etall (|200lD estima te 0.14t°;^ 
with e rrors at 95% confidence, for SDSS. iSzalav et alJ 
(2001) also give an estimate of as of 0.92 ± 0.06, but 
the issue of galaxy biasing is folded into this determina- 
tion. Although there are indications from the 2dF survey 
that biasing for the relevant galaxies is nearl y unity from 
redshift-space distortions (b = 1.04 ± 0.11. IVerde et alJ 
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Fig. 3. — The prior probability used for the shape parameter 
r e ff is shown as solid gr een. (This is to be contrasted with the 
more skewed one used in Lan ge et alJ (2001), etc., shown as light 
dashed-dotted green.) The T-prior was designed to encompass the 
range indicated by the APM data, vintage 1992, but estimates 
from the 2dFRS and SDSS shown below are quite compatible. The 
CMB results for F g q are shown with various choices for priors for 
comp arison. Sol i d blue is for 'all-data' as of June 2002, as described 
in Sievers et al. (2003), and magenta is for the data as of March 
2003, as described in Bond et al. (2003). T e ff includes corrections 
for h, and the tilt. The related values of Qmh are shown at the 
bottom to show the effect of these corrections. 

2002), and the result is therefore compatible with the re- 
sults shown, such results are not as directly applicable 
as those from clusters and from lensing which directly 
relate to matter density power spectrum amplitudes. 

The LSS prior used here and in Paper V involves 
a combination of constraints on the amplitude param- 
eter crg^^n 56 an d on the shape parameter r W. It 

differs slig h tly fr om t hat used in lLange et all ([200 It) . 
Jaffe et all J2001) and iNetterfield et all (|2002D . We use 
ct8^^ 56 =0.471q"o2'1o.08' distributed as a Gaussian (first 
error) smeared by a uniform (tophat) distribution (sec- 
ond error). We constrain the shape of the power spec- 
trum via r eff = T + (n s - l)/2=0.2llo o3'+o og, where 

T * ^k-^^ 1 ^ 1 ^" ' 06 ]. In earlier work a 
central value of as^t^ 56 = 0.55 was used together with 
Teff = 0-22^o!o4 io'.o?- The T prior is compared with 
the current version in Fig. |3J The change does not affect 
the parameter values obtained. 

The original as prior choice was motivated by fits to the 
cluster temperature distribution, as was the decision to 
lower the central value by 15% made here. Our philoso- 
phy is to make the distribution broad enough so that rea- 
sonable uncertainties are allowed for in the prior. For ex- 
ample, there are many models that do not fit the shape as 
well as the amplitude of the cluster distribution function. 
Thus the best as for a given model depends upon the 
temperature range chosen for the fit, and other physics 
might be involved. Especially with the reduced model 
spaces often considered, the formal statistical errors can 
look spectacularly good, but systematic issues undoubt- 
edly dominate. Curiously, with the 15% drop, it appears 
that the as prior chosen could have been designed for the 
(weak) weak lensing results. We test sensitivity by show- 
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Fig. 4. — One and two sigma contours of 2D projected likelihood 
functions show how various cosmological parameters correlate with 
(78- For this case, 'all-data' from the June 2002 compilation was 
used: DASI+MAXIMA+BOOMERANG+VSA plus the CBI mo- 
saic data for the odd AL=140 binning. The Qj~ — as panel shows 
the weak-h prior (magenta) and LSS+weak-/i prior (solid black). 
For the other 3 panels the flat constraint was added to these two 
priors as well. Note the positive correlation with n s , but little 
correlation in the other variables. 

ing results for a 'LSS (low-as)' prior which has a further 
drop of 15%: c r 8^m 56= 0-40lQ"o2'lo.08' This also accom- 
modates some of the recent lower as estimates from clus- 
ter abundances, many of which use the X-ray luminosity 
as an indicator of mass, calibrated by observations. (We 
note that the re-anal yzed V irmos-Descart Survey result 
(van Waerbeke etMlUM) gives a 8 ^ 56 =0.43 ± 0.04, 
similar to the value obtained from the CMB-only data 
when WMAP is included.) 

We will discuss the simulations in detail below, but 
we note here that the cases we have run simulations on 
have a 8 = 0.9 and 1.0, and r eff = 0.18 and 0.21. The 
detailed baryonic dependence of the transfer function was 
included in the 0.18 case. The parameter choices were 
cj b = 0.022 and 0.0245, h = 0.7, tt m = 0.3 and 0.37. A 
best-fit model to the data for many prior choices has uj b = 
0.0225, u c = 0.12, n s = 0.95, r c = 0.1, and Q A = 0.7 and 
with h = 0.69, Q b = 0.047, ^ cdm = 0.253 with a 13.7 Gyr 
age; the normalization is as = 0.84, and T = 0.17 with 
r eff = 0.15. 

In general, CMB data provide weak constraints on the 
normalization of the matter fluctuations due to the added 
effects of the spectral tilt n s and the optical depth param- 
eter rc on the overall amplitude of the CMB power spec- 
trum. Allowing for tensor modes in the power spectrum 
introduces further degeneracies. As shown in Paper V, 
when fitting for a given parameter we adopt the conser- 
vative, and more computationally intensive, approach of 
marginalizing over all other variables considered. 

In Table we summarize the constraints on the LSS 
parameters from the CMB data used in Paper V and a 
combination of priors. The mean and quoted (la) errors 
correspond to the 50%, 16% and 84% integrals of the 
marginalized likelihoods respectively. 

In Figs. E] and [HI we show the correlation of as with 
other parameters considered. Correlations with rc are 
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TABLE 1 

Simulation parameters 



Code 


Size (Mpc) 


Resolution 


^8 


n b h 2 


^A 




h 


n s 


r e ff 


SPH 


200 


256 3 


0.9 


0.0200 


0.70 


0.30 


0.70 


1.0 


0.18 


SPH 


200 


256 3 


1.0 


0.0200 


0.70 


0.30 


0.70 


1.0 


0.18 


SPH 


400 


512 3 


0.9 


0.0220 


0.70 


0.30 


0.70 


1.0 


0.18 


MMH 


143 


512 3 


1.0 


0.0245 


0.63 


0.37 


0.70 


1.0 


0.21 



TABLE 2 

Amplitude and Shape Parameters for LSS from the CMB: 'all-data', June 

2002 Compilation 



Priors 



^8 



^8 



^ 56 



- eff 



wk-/i 

wk-/i+LSS 
wk-/i+LSS (low-<78 ) 
wk-/i+SN 



0.88; 
0.83; 

o.8i; 

0.85; 



0.61016 
0.52° 
0.45° 
0.45g 



n 28 - 07 

Q-QT 
0.230 

0.21° 

0.2lg 



0.230-06 
0.20| 
0.17| 
0.17! 



0.210-05 

0.18! 
0.17| 
0.17! 



' 



03 



flat+wk-/i 


0.892 


.10 

■hi 


0.560 


23 


0.26°; 


flat+wk-/i+LSS 


0.87° 


0.500 


8? 
hh 
?8 

11 


0.23°; 


flat+wk-/i+LSS(low-cr 8 ) 


0.84° 


:?! 

.09 


0.43° 


0.20°; 


flat+wk-h+SN 


0.85g 


0.43g 


0.20°; 



0.232 
0.20S 
0.172 
0.17° 



0.202 
0.17° 
0.15° 
0.15g 



flat+HST-Ji 


0.860 


11 

hi 
?f 


0.44° 


13 

41 


0.21S 


flat+HST-/i+LSS 


0.880 


0.47° 


0.222 


flat+HST-/i+SN 


0.840 


0.43° 


81 


0.202 


flat+HST-/i+LSS+SN 


0.87° 


81 


0.45° 


8? 


0.21° 


flat+HST-/i+LSS(low-cr 8 )+SN 


0.85° 


81 

07 


0.43° 


8§ 

05 


0.20° 



0.18° 


05 

m 


0.16° 


0.19° 


0.17° 


0.17° 


0.15° 


0.18° 


81 


0.16° 


0.17° 


81 

02 


0.15° 



Note. — Amplitude and shape parameter estimates from various datasets, for various 
prior probability choices and using 'all-data' for the CMB available for Paper V, and as 
described there, in particular DMR+DASI+MAXIMA + BOOMERANG (Ruhl et al. 2003 
cut), VSA and CBI mosaic data for the odd A^=140 binning, as described in Paper III. 
In the first four rows, the weak prior in h (0.45 < h < 0.90) is imposed (including further 
weak constraints on cosmological age, £0 > 10 Gyr, and matter density, O m > 0.1). The 
sequence shows what happens when priors for LSS (with agQ^ 56 centered about 0.47), for 
LSS (low-as) (with the ag distribution shifted downward by 15%, centered about 0.40), and 
for SN are imposed. While the first four rows allow Qtot to be free, the next four have fltot 
pegged to unity, a number strongly suggested by the CMB data. The final 4 rows show 
the 'strong-ft,' prior, a Gaussian centered on h — 0.71 with dispersion ±0.07, obtained for 
the Hubble key project. The to > 10 Gyr and Q m > 0.1 constraints are also imposed, but 
they have no impact. Central values and la limits for the 7 database parameters that form 
our fiducial minimal-inflation model set are found from the 16%, 50% and 84% integrals 
of the marginalized likelihood for ag and agft,^ 56 . For T eq = Q m h, T and r e ff , the values 
are means and variances of the variables calculated over the full probability distribution. 
Shifting the center of the agQ^ 56 prior downward by 15% to accommodate more of the 
low cluster results has little effect on the flat+HST-ft,+LSS+SN result. 



nontrivial, and allowing the space to be open in rc does 
tend to give higher marginalized The CMB data as 
of June 2002 did not give an indication of a rc detection, 
and we have imposed no physical prior on those results. 

The addition of the first-year WMAP data and other 
data that has appeared since sharpens the determina- 
tion of the parameters as, <T8^^ 56 and r e ff. This 'March 
2003' comp i lation of CMB anisotropy data, described in 
iBond et al.l ((2003), is treated in the same way as the June 
2002 data is in Paper V, with the same Q-database. Re- 
sults are shown in Table [31 and in the reduced errors in 
Figs, m [21 and [31 For the March 2003 data, a prior was in- 
cluded on rc to accommodate the rc = 0.16 ±0.04 detec - 
tion reported by the WMAP team (jKoerut et al.1 12003^. 
The form we have adopted has a top-hat spread con- 
volved with a Gaussian distribution, as for our crg^^ 56 



and r e ff priors: rc = 0-16^q"o4+^06- Similar results are 
obtained using the Markov Chain Monte Carlo approach, 
which explicitly includes t he rc-detection using the TE 
data of lKogut etap (|2003ft . 

Best-fit models may have lower as values than the 
marginalized values. Typical best-fit models for the 
June 2002 dataset are: (D cdm = 0.47, ^ A = 0.50, tt b = 
0.079, h = 0.51, n s = 0.90, r c = and a 8 = 0.72) for 
the weak prior case and (£l c dm = 0.44, = 0.50, ftb = 
0.063, h = 0.57, n s = 0.90, r c = and a 8 = 0.81) for 
flat+weak and flat+weak+LSS priors. These a 8 values 
are lower because the best-fits selected the rc = peak 
in the likelihood whereas integration over the relatively 
broad rc likelihood allows for the inclusion of high-ag 
models. If the HST-/i or SN prior is included along with 
LSS, the best-fit model is a more conventional ACDM 
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TABLE 3 

Amplitude and Shape Parameters for LSS from the CMB: 
March 2003 Compilation 



ALL-DATA . 
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Note. — Similar to Table. EJ but for an extension of 'all-data' to include WMAP, 
ACBAR, extended VSA ( ( Grainge et al. I 12003^ and CBI results, as well as DASI, 
Boomerang and Maxima, analyzed in the same way. A further (weak) prior on rc was 
applied to accommodate the WMAP detection from the cross-correlation of temperature 
and polarization, but this makes little difference to the results. For example, without this 
extra prior, flat+wk-/i shifts to 0.81°;°5 and 0.4lg;??. As for Table^l the downshifted LSS 
prior has little effect on flat+HST-/i+LSS+SN. 




Fig. 5. — One and two sigma contours of 2D project e d like lihood 
functions for uj h and tq for 'all-data' of Sievers et al. (2003). The 
flat+weak-ft, prior (magenta) and flat +LSS+ weak- h prior (solid 
black) cases are shown. The significant ctq-tc correlation is ev- 
ident, which results in a higher as for higher tq- One could im- 
pose a stronger prior than rc < 0.7 (as is done here) based on 
astrophysical arguments. This is fraught with uncertainty since it 
involves the first objects collapsing on small scales in the universe 
and their efficiency in generating stars that produce ionizing radi- 
ation. However, rc has apparently bee n detected by WMAP at 
the 0.16 ± 0.04 level (JKogut et al.l l2003h and results with a prior 
encompassing this detection on the March 2003 compilation of the 
data are given in Table El 

one, with rc = 0.1: (^ c dm = 0.25, = 0.70,^6 = 
0.047, h = 0.69, n s = 0.95, r c = 0.1 and a 8 = 0.84). For 
the March 2003 data with the tq prior, this rc = 0.1 
model also provides a best-fit for the weak+flat and 
weak+flat+LSS cases. 



Comparison of our results with independent estimates 
of erg from CMB data is complicated by the different 
choices of parameter margin alization and by different 
treatments of the CMB data. lLahav et alJ (|2002l) carried 
out a joint CMB-2dFRS analysis of cosmological param- 
eters. Their result is erg = 0.73 ±0.05 after marginalizing 
over some variables but keeping rc fixed at and n s fixed 
at unity. This value is therefore biased low with respect 
to non-zero rc models. They also show how lett i ng rc 
vary gives higher values for ag. Melchi orri fc Silkl {2002) 
quote smaller values for in their analysis, even when 
rc was allowed to vary. 

The values we obtain for the March 2003 dataset are 
in go od agreement with those obtained by the WMAP 
team (Sne rgel et alll200l . However we caution that the 
projected distribution of as from the CMB depends upon 
the other parameters. For the seven-parameter inflation- 
motivated models considered here, there is a strong corre- 
lation of erg with the Thomson scattering depth rc and 
the primordial spectral index n s . If, in addition, n s is 
allowed a logarithmic correction with wave number, in- 
troducing another parameter highly correlated with 
the distribution extends to higher as ( Bond et al. 2003). 
Allowing for a gravity wave contribution extends the dis- 
tribution to lower ag. These caveats about the effect of 
adding extra parameters in the CMB analysis together 
with the scatter in the estimates of as from clusters and 
weak lensing shown in Figs. H an d El reflect the current 
uncertainty as to whether a low or high as will emerge. 

3. THE ANGULAR POWER SPECTRUM OF THE 
SUNYAEV-ZELDOVICH EFFECT 

The SZ effect is a signature of the scattering of CMB 
photons off hot electrons. The effect can be described in 
terms of the fractional energy gain per scatter along the 
line of sight. By multiplying by the number density of 
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electrons and integrating along the line of sight we can 
derive the induced temperature change: 



AT SZ 

^CMB 



-2y^ K {x) 



-2a T 



kB(T e — Tcmb) 



m e c z 



d\ ^k(x) 
(1) 



where y is the Compton ^/-parameter, gt is the Thomson 
cross-section, n e is the electron number density, T e is the 
electron temperature, Tcmb is the CMB temperature, 
ks is Boltzmann's constant and m e c 2 is the electron's 
rest mass energy. Here i/jk(x), % = hv/ksTcMB is a 
frequency-dependent function which is unity at Rayleigh- 
Jeans wavelengths and is 0.975 at the 30 GHz frequencies 
probed by the CBI. 

With the great increase in experimental sensitivity, 
measurement of the SZE in known clusters has become 
almo st routine (see, e.g. , [Birkinshaw. Gull & Hard ebeckl 
" 19841: iCarlstrom et all 1199ft iHolzanfel et all 1 1997: 



Mason. Mvers. fc Readheadl 120011: lUdomorasert et all 



2004: iGrainge et al. 112002ft . Coupled with X-ray obser- 
vations of the cluster these measurements yield indepen- 
dent constraints on the value of the Hubble constant. 
They also provide a direct measurement of the amount of 
baryon mass in the cluster gas. The SZ effect is expected 
to contribute significantly to the CMB power spectrum 
at scales i > 2000, with a crossover point between the 
primary CMB and SZ signature occuring somewhere be- 
tween i ^ 2000 and i ~ 3000. Surveys observing at these 
scales will therefore require accurate component separa- 
tion to reconstruct the primary CMB spectrum and much 
work has focused on this issue in recent years. Many pro- 
posed experiments will adopt multi-frequency observing 
strategies in order to separate the different components 
by their different spectral dependence. These techniques 
cannot be applied to the narrow frequency band of the 
CBI observations. However we can attempt to address 
the question of whether the SZ effect could provide a 
contribution with the required amplitude to explain the 
observed excess. 

We use the output of two separate hydrodynamical 
simulation algorithms in an attempt to predict the level 
at which the SZ effect will contribute to the CBI deep 
field observations. We also relate our numerical results 
to analytical models of the SZ power spectrum based on 
the halo model. 

3.1. Hydrodynamical Simulations of the 
Sunyaev-Zeldovich Effect 

The two codes that we have employed both provide 
high-resolution dark matter and hydrodynamical sim- 
ulations of large scale cosmic structure which we use 
to generate simulated wide-field SZ maps. The simu- 
lation algorithms an d processing techniques were devel- 
oped independ ently (iPenlll 998 -aUBond et, al .111991 2002: 
nnl2 



Wadslev. Stadel fc Ouinnl 2003) and are based on two dif- 
ferent numerical schemes for solving the self-gravitating, 
hydrodynamical equations of motion. High-resolution 
gasdynamical simulations of large enough volumes are 
still beyond current technological capabilities. Follow- 
ing an approach taken previously b y other authors (e.g., 
iSpringel. White fc HernauTstll2001l) . we create a pseudo- 
realization of the cosmic structure up to high redshifts 
from a single, medium-sized, high-resolution simulation 
by stacking randomly translated and rotated (or flipped) 



copies of the (evolving) periodic volume. 

One set of SZ simulations was obtained using the 
Gasoline code, an efficient implementation of the (La- 
grangian) smooth particle hydrod ynamics (SPH) method 
' (fWadslev. Stadel fc Ouinll200l . This tree+SPH code 
uses a pure tree-based gravity solver and has spatially- 
adaptive time stepping. It has been parallelized on many 
architectures (MPI in this case) and shows excellent scal- 
ability. The results presented here are based on the anal- 
ysis of three high resolution ACDM simulations: two 200 
Mpc box computations, with 256 3 dark matter plus 256 3 
gas particles, and one 400 Mpc box computation, with 
512 3 dark and 512 3 gas particles, a very large number 
for SPH simulations. The calculations were adiabatic, in 
the sense that only shocks could inject entropy into the 
medium. Despite the different sizes, all three simulations 
were therefore run with the same mass resolution. All 
simulations were run with a gravitational softening of 50 
kpc (physical). The scale probed by the SPH smoothing 
kernel was not allowed to become smaller than the grav- 
itational softening in the 200 Mpc runs, but it was not 
limited in the 400 Mpc run, which attained gas resolu- 
tion scales as low as 5 kpc in the highest density environ- 
ments, though increasing considerably at lower density. 
(64 neighbors were required to be within the smoothing 
kernel in the 400 Mpc run.) The 512 3 simulation was 
performed on a large- memory, 114 (667 MHz) processor 
SHARCNET COMPAQ SC cluster at McMaster Uni- 
versity. It required about 80 GB of memory and took 
~ 40 days of wall time to run. 

Another set of SZ maps was obtained using a 512 3 
run of the M oving Mesh Hydrodynamics (MMH) code of 
iPenl (|1998al) . This code features a full curvilinear Total- 
Variation-Diminishing (TVD) hydro code with a curvi- 
linear particle mesh (PM) N-body code on a moving co- 
ordinate system. We follow the mass field such that the 
mass per unit grid cell remains approximately constant. 
This gives all the dynamic range advantages of SPH sim- 
ulations combined with the speed and high resolution 
of grid algorithms. The box size was 143 Mpc and the 
smallest grid spacing was 57 kpc (comoving). This 512 3 
simulation used 30 GB memory and took about three 
weeks (~ 1500 steps) on a 32 processor shared memory 
Alpha GS320 using Open MP parallelization directives. 
The calculations were adiabatic as well. 

The simulation boxes yield a number of projections of 
the gas distributions in random orientations and direc- 
tions. The projections are then stacked to create a red- 
shift range appropriate for the SZ simulation. In the case 
of the MMH simulation the angular size of the simulated 
box above a redshift of z = 1.6 is smaller than 2 degrees, 
the required angular size of the SZ simulations. To in- 
clude the contributions from higher redshifts we tiled two 
copies of the periodic box to cover the increased angles. 
This only affects the smallest scales. All the maps were 
obtained with a low redshift cutoff at z = 0.2 to mini- 
mize the contribution from the closest clusters in the pro- 
jections. The computational costs of running such large 
hydrodynamical simulations prevented us from obtaining 
targeted simulations with identical cosmological param- 
eters from both algorithms. Below we account for the 
parameter variations in the models run when comparing 
the two simulations and we see a remarkable consistency 
between the results. The parameter variations are, in 
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Fig. 6. — SZ power spectra for various simulations. All SZ 
spectra have been scaled to a common Q^h = 0.0314, the best- 
fit to the CMB data. The upper panel shows the SZ spectra 
for the cr 8 = 1.0 SPH 256 3 200 Mpc (cyan) and MMH 512 3 140 
Mpc (red) simulations. We also plot the power spectrum from 
the a 8 = 0.9 SPH 256 3 200 Mpc run (yellow points) scaled using 
the relation of equation The target 2-sigma region suggested 
by the deep data is shown as a green-hatched box. An optimal 
power spectrum combining all CBI data with the BOOMERANG, 
COBE-DMR, DASI and MAXIMA data is shown as blue trian- 
gles. The solid black curve is a best-fit model to the data out to 
t = 2000. The middle panel compares the a 8 = 0.9 512 3 400 Mpc 
(cyan) and 256 3 200 Mpc (magenta) SPH simulations. The yel- 
low points are the spectra derived from the ACDM simulation of 
SpringeL White XHernauisti (l200lL also scaled to Q b h = 0.0314. 
The bottom panel compares the MMH results (dot-dash) with 
the analytic halo model results described in § |3| for mass cuts of 
O.OOlMs (solid, red), 0.01M 8 (dotted, blue) and 0.1M 8 (dashed, 
magenta) described in the text. There is little sensitivity to a 
change in the lattice size over the Grange shown. We note that 
physical effects, such as early entropy injection, may change the 
effective mass cut. 

fact, an advantage for us since they provide us with a 
sampling of the SZE amplitude for processing through 
the CBI pipeline. We summarize the parameters used in 
the simulations in Tabled 

The two 200 Mpc SPH simulations were used to gen- 
erate 20 2x2 degree maps each while the MMH simula- 
tion yielded 40 separate maps. A detailed analysis of the 
MMH map s including SZ statisti c s and non-Gaussianity 
is given by IZhang. Pen fc Wand ()2002t) . A preliminary 
analysis of simulated SZ maps base d on one of the 256 3 
SPH simulations was presented bv iBond et alJ (|2002fi . 
and a more detailed analysis of the 512 3 results will be 
presented elsewhere. 

3.2. Hydrodynamical Simulation Results 

In Fig. El we show the results for the os = 1.0, 143 
Mpc MMH simulation, the cr 8 = 1.0 and 0.9, 200 Mpc 
SPH simulation and a 8 = 1.0, 400 Mpc SPH simula- 
tion. The curves show averages for the Ci spectra from 
the 40 and 20 MMH and SPH maps respectively. The 



shaded regions show the full excursion of the power spec- 
tra in each set of maps. This gives a partial indication 
of the scatter induced in the power due to the small ar- 
eas being considered. We compare these with a best-fit 
model to the BOOMERANG, CBI, COBE-DMR, DASI 
and MAXIMA data and an optimal combination of the 
mosaic and deep CBI, BOOMERANG, COBK-DMR, 
DASI and MAXIMA bandpowers (Paper V). The com- 
bined spectrum is designed for optimal coverage with 
variable bandwidths over the range of scales considered. 
At scales above £ ~ 1000 the optimal spectrum is domi- 
nated by the contributions of the CBI mosaic and deep 
results. 

To compare the cosmological simulations we have to 
account for the difference in parameters the simula- 
tions were run with. The dominant effect for the SZ 
Ct spectrum is variation in os and fl^h, scaling as a| 
and (Qbh) 2 . The top panel of Fig. demonstrates that 
Ct oc crj does indeed bring the 200 Mpc SPH run at 
erg = 0.9 into essentially perfect alignment with the 
erg = 1 run (yellow/square points). (Apart from the am- 
plitude, the initial conditions in the simulations were oth- 
erwise the same). The exact factor n in a scaling of form 
Ct (Q b h) 2 (Tg depe nds upon the model in ques tion, but 
6 < n < 9 is typical (IZhang. Pen fc WandEool. Differ- 
ences in Q m and shap e of the power spec t rum a lso have 
an influence on C^, as iKomatsu fc Sell aid (|2002ft discuss 
in more detail. 

For the two separate <7s cases employed in the runs we 
rescale the spectra to a nominal Q^h = 0.0314, a value 
suggested by the CMB data. This enables us to compare 
the MMH result with the SPH results for the os = 1 runs 
(top panel Fig.EJ). We see that the two codes give similar 
amplitudes over the range considered, and, in particular, 
these are very similar over the scales of interest for the 
CBI deep field result (green/hatched box). At larger 
angular scales, the MMH simulation shows a somewhat 
higher amplitude than the SPH simulation. This may be 
in part due to the presence of a large rare cluster in the 
MMH volume simulation. The Poisson noise contribu- 
tion of the large cluster is aggravated by the resampling 
technique adopted in the SZ map making stage described 
above. Variable redshift cutoffs do indeed confirm the 
dominant effect of this single cluster on the low-^ tail of 
the spectrum. 

At smaller angular scales the two results start to di- 
verge. The differences in spectral shape and cosmolog- 
ical parameters should not account for this. It may be 
attributable to the different techniques used to limit the 
resolution achievable, but this needs more investigation. 

In the middle panel we have plotted the re- 
sults from the 400 Mpc and 200Mpc SPH simula- 
tions which were both run with os = 0.9. We 
also compare this run to SZ power spectra de- 
rived from the 200 Mpc box SPH simulations of 
iSpringel. White fc Hernauistl (120011). The cosmologica l 
parameters of the SpringeL White & Hernquist (2001) 
simulations were the same as those used in our SPH runs 
except for a value of Vt^h 2 — 0.018. We also rescale this 
spectrum to the nominal Q^h = 0.0314. The agreement is 
remarkable. We can contra st this level of agreement with 
the si tuation described in iSpringel. White fc HernouTsU 
(2001) where it appeared that different codes were giv- 
ing quite different results. Nonetheless we plan a further 
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exploration of the effects of lattice size variations and 
other numerical parameters to compare more exactly the 
MMH and SPH codes. 

Based on these results we can calibrate the expected 
power to compare with the wide band result of Paper II: 



SPH: Cf z ~ 170 fiK 2 



(0.0314) 2 



MMH: Cf z ~ 210 ^K 2 ~ 7 



(2) 



(0.0314) 2 8 * 

For a 8 = 1 and Q^h = 0.0314, these amplitudes are ~ 0.4 
relative to the noise level for the CBI joint analysis of the 
three deep fields. 

3.3. Analytic Modeling of the Sunyaev-Zeldovich Power 

Spectrum 

The general analytic framework comes under the 
generic name of 'halo models', in which simple pa- 
rameterized gas profiles within clusters and groups 
are constructed, appropriately scaled according to the 
masses. The gas profile in a halo and halo mass- 
temperature relation determines the SZ effect of the 
halo. The abundance of the halos as a function of mass 
and r edshift is dete rmined by the iPress fc Schechteil 
(|1974fl formula, the iBond fc Mverl (|1996l) peak-patch 
formula, or formulas derived from fits to TV-body simula- 
tions. Clustering of the halos is included through simple 
linear biasing models. The halo model approach has 
been applied to th e SZ effect by many a uthors over time 
(e.g.. iBonl H983 iCole fc Kaiserl Il983 iMakino fc Sut< ' 
1993; IBond fc Mversl Il996t lAtrio-Barandela fc Mucketl 
1999: iKomatsu fc Kitayamal 119991: iCooravl 1 2000: 
Mol nar fc Birkinshawl 120001: LSeliak. Burwell fc Perl 
2001). There has still not been an adequate calibra- 
tion of the relation between the physical parameters 
describing the gas distribution in the halos and the 
results of gas dynamical simulations, and usually the 
parameters that have been adopted are those derived 
from observations of clusters at low redshift. We show 
here that we can get reasonable fits to our simulation 
results by choosing suitable gas profile parameters. 
However we caution that the gas profile is a function 
of halo mass and redshift, so globally fitting to our 
simulated SZ power spectra may not deliver reliable 
parameters for the analytic model. For now, we take 
our globally fit models to show sensitivity to parameter 
variations. 

In the bottom panel of Fig. El a few analytic models are 
shown to compare with the <r 8 = 1 MMH simulation. The 
Press- Schechter distribution for the halo comoving num- 
ber density n as a function of halo mass M and redshift 
z is Mdn/dM oc (2/tt) 1 / 2 (po/M) \d\ncj/d\nM\ ve'^l 2 , 
where v = (S c /a). Here po is the present mean matter 
density of the universe. <r(M, z) is the linear theory rms 
density fluctuation in a sphere containing mass M at 
redshift z. It is calculated using the input linear density 
power spectrum P(k) of our simulations, which is trun- 
cated at /ci ower = 2tt/L and k upper = ^ttN/L due to the 
finite box size L and resolution N of simulations, respec- 
tively. In this formula, we have taken S c = 1.686 as the 
linearly extrapolated spherical overdensity at which an 
object virializes. Although this is strictly valid only for 
Q m = 1 cosmologies, when Q m decreases from 1 to 0.3, 



S c only decreases from 1.686 to - 1.675 ()Eke et al.ll 996). 
We omit this dependence of 5 C on Q m , and therefore 
redshift, for the Q m = 0.37 ACDM we are trying to fit 
here. We also truncate the mass integral at a specified 
lower mass limit Mi ower , which we choose as a free pa- 
rameter. In these analytic estimates, clustering is usually 
treated with simple linear models in which the long wave- 
length cluster distribution is amplified over the underly- 
ing mass distribution by a mass-dependent biasing factor, 
leading to a power spectrum b(M\)b(M2)P (k) amplified 
over the underlying dark matter spectrum P(k). This 
turns out to be a small effect for the SZ power spectra. 

To describe the gas profile, we need the pressure pro- 
file. We take this to be that for an isothermal distribu- 
tion with a baryon density profile given by a f3 model 
with (3 = 2/3. For the gas te mperature, w e adopt the 
virial theorem relation given in QE3 (11998b! ) : M/M 8 = 
(T/T 8 ) 3 / 2 , where M 8 is the mean mass contained in an 

Sfc^Mpc sphere and T 8 = 4.9^ /3 ft(z)°- 283 (l + z) keV. 
Q(z) is the fraction of matter density with respect to 
the critical density at redshift z. This relation was ob- 
tained by comparing the gas temperature distribution 
in simulations with the halo mass function described 
above. For the electron number density profile we adopt 

n e = n e o [l + r2 / r core] •> scaled by the central density 
n e o and core radius r core . Since our simulations show 
that the pressure profile falls off more rapidly than the 
product of temperature and density given here would in- 
dicate, we truncate the pressure at a fraction /r v i r of 
the 'virialized radius' r v i r . The virial radius is defined 
as the radius of a sphere with mass M and mean den- 
sity A c (z)p(z), where p(z) is the mean matter density at 
redshift z and A c is given bv lEke et alJ ((199 6). If we as- 
sume that gas accounts for a fraction VLb/ ^™ of the halo 
mass, the baryon content fixes one of the 3 parameters, 
n e o- We treat / and r v i r /r core as free parameters to be 
fit to the SZ power spectrum. 

The bottom panel of Fig. El shows typical fits we can 
obtain with the above model using r v i r /r core = 4.8 and 
/ = 0.9. Given the number of free parameters, many 
combinations of values can yield reasonable fits. In or- 
der to investigate the resolution effects in our simulations 
we fix r v i r /r c and / and vary the lower mass cutoff Mi ower 
and the k range cutoff in the analytic model. These pa- 
rameters are related to the resolution limitations of our 
simulations. For the MMH simulation, we find the effect 
of the k range cutoff is negligible, but that Mi ower has a 
larger effect. Since we need more than 100 gas and dark 
matter particles to resolve a halo, we choose Mi ower as 
the mass of 100 gas particles and 100 dark matter par- 
ticles. This corresponds to 0.0015M 8 for the 512 3 sim- 
ulation. We see substantial deviations developing when 
Slower > O.OlMg. There are as well other uncertainties 
not included in our hydro simulations, which are adia- 
batic. For example, the outflow of gas from groups could 
lead to an effective mass cutoff that is physical rather 
than resolution dependent. The concentration parame- 
ter for gas, r v i r /r cor e could be a function of mass and of 
redshift, which would also change the results. 



4. 



SIMULATING CBI OBSERVATIONS OF THE 
SUNYAEV-ZELDOVICH EFFECT 
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The estimation of bandpowers used in the analysis of 
the CBI observation is based on a maximum likelihood 
technique which assumes the signal to be a Gaussian ran- 
dom field (Paper IV). This assumption appears to be 
justified in the case of primary CMB anisotropics where 
no clear evidence of non-Gaussianity has been found. 
However, the contributions from various foregrounds are 
non- Gaussian. Although the observations are noise dom- 
inated at the scales of interest for the CBI deep field ex- 
cess and therefore expected to be predominantly Gaus- 
sian, it is important to test the effect of a small non- 
Gaussian signal such as a SZ foreground on the band- 
power estimation procedure. 

We used the simulated SZ maps to test the effect of a 
SZ foreground on the CBI analysis pipeline. To do this 
we constructed detailed mock data sets using the CBI 
simulation tools (Paper IV). We took real observations 
for particular dish configurations and pointings and re- 
placed the observed visibilities by realizations of the ex- 
pected signal and instrumental noise. This resulted in 
simulated observations with uv coverage identical with 
the actual observations. The simulations can also include 
foreground templates as a distribution of point sources 
or maps of the SZE. 

Each simulated 2x2 degree SZ map is used as a fore- 
ground and added to a primary CMB background to gen- 
erate mock observations of the 08 h deep field. One point 
to note is that the CBI observations are actually differ- 
enced to minimize any ground pick-up (Paper II). This 
involves subtracting the signal from two fields (lead and 
trail) separated by 8 minutes of right ascension on the 
sky. In adding a SZ foreground to the simulated obser- 
vations two separate SZ maps were used as lead and trail 
fields. For this exercise, we created 20 simulated obser- 
vations using the 40 MMH maps and 10 for the 20 SPH 
maps. 

The mock data sets were then processed through our 
power spectrum estimation pipeline described in Pa- 
per IV. The maximum likelihood calculation of the power 
spectrum used template correlations to project out the 
expected contribution from point source foregrounds. 
Since point source templates were projected out of the 
data or were constrained at a known amplitude (as in the 
case of the residual unresolved background) , the power in 
the excess was assigned to the CMB when fitting for the 
bandpowers since no other correlations were included in 
the quadratic power spectrum estimator (except for in- 
strumental noise). We reproduced this situation by sim- 
ulating observations of primary CMB realizations with 
SZ foregrounds, but only allowing for noise and CMB 
contributions to the total correlation matrix in the like- 
lihood 

C = C N + C CMB , (3) 

where the correlations due to the CMB are expressed as 
a function of the bandpowers Cb with 



iCMB 



dC 



CMB 



dC f 



(4) 



The SZ maps were produced from simulations of dif- 
ferent cosmologies and we have shown how the spectra 
can be scaled to fiducial Q^h and as values for com- 
parison. However when simulating the measurement of 




Fig. 7. — Simulated observations of CMB plus SZ fields. The 
result of bandpower estimation on simulated observations of the 
08 h deep field containing noise, realizations of a fiducial ACDM 
model shown in black (solid line) and SZ foregrounds from the SPH 
and MMH erg = 1.0 simulations. The blue (dashed) curve shows 
the average power spectrum of the input SPH maps and red (dash- 
dotted) is the average power spectrum of the input MMH maps. 
The red (triangle) points are the average bandpowers obtained from 
the analysis of 20 separate observations using the MMH maps. The 
blue (square) points are the average bandpower obtained from 10 
SPH maps. The bandpower estimation pipeline recovers the correct 
power in the SZ dominated region t > 2000. The green (shaded) 
rectangle shows the 95% confidence region for the high-^ CBI deep 
field result. Both the MMH and SPH codes show power consistent 
with the confidence region. 



bandpowers we chose to use unsealed contributions in or- 
der to test the different amplitude regimes given by the 
different values of baryon density. As a primary CMB 
signal we used random realizations of a single ACDM 
model (Q m 0.30, Q b 0.04, Q A 0.70, h 0.68 
and n s = 0.975). For the SZ contributions we restricted 
ourselves to the as = 1.0 simulations. 

Results are shown in Fig. The bandpowers are av- 
eraged over the 20 and 10 independent realizations from 
the two as = 1.0 simulations. The triangles (red) are for 
the MMH simulations and the SPH results are shown as 
squares (blue) points. The model used for the primary 
CMB contribution is also shown together with the 2a 
confidence region for the deep excess. The errors shown 
are obtained from the variance of the measured bandpow- 
ers. The averages appear to recover the input power ac- 
curately in both regimes where either the primary CMB 
signal or SZ signal dominate. 

The interpretation of the errors is of course com- 
plicated in the SZ-dominated r egime by the non- 
Gaus sianity of the SZ signal (eg. IZhang. Pen fc W ang 
2002). The small area considered in the simulations re- 
sults in significant sample variance effects in the measure- 
ments of the SZ power. The scatter in the observed band- 
powers over the different realizations is also shown in 
Fig. [3 The non-Gaussian scatter is significant, although 
this effect would decrease for the larger area probed in 
the joint three-field analysis and because the high-^ er- 
rors are dominated by the noise. It is also important to 
note that if the SZE is a significant source of the power in 
the observations, the sample variance component of the 
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errors derived in the optimal spectrum estimation would 
be biased. 

It is also interesting to note the effect of the \ow-£ 
contribution in the MMH simulations. There is a siz- 
able contribution to the bandpowers below t ~ 2000. 
Although this effect may be due to the presence of one 
rather large cluster in the simulation, as discussed above, 
it raises the question of whether a SZ signal with high 
enough amplitude to explain the excess may already be 
constrained by its Vm-t contributions. 

5. WIENER-FILTERING THE CBI DEEP FIELD 
OBSERVATIONS 

The CBI does not have sufficient frequency coverage to 
strongly distinguish different signals by their frequency 
dependence (e.g., CMB, SZ, and Galactic foregrounds). 
Spectral separation of the different components is there- 
fore infeasible. Although the sensitivity of the CBI re- 
sults presented here is not sufficient to conclusively iden- 
tify individual features in the maps, we present examples 
and simulations of a Wiener- filtering technique described 
in Paper IV. 

The total mode to mode correlations in the data are 
separated into four components in a typical analysis of 
the deep field visibilities with the weight matrix W as 
the sum 

W = (C N + C CMB + C SRC + C res )-\ (5) 

with noise (N), CMB, known source (SRC) and resid- 
ual source (res) components, respectively. (In papers II 
and III, the known source contribution is split into two 
components, for sources with and without measured flux 
densities.) 

In the limit where all the components can be consid- 
ered Gaussian random fields we can define the prob- 
ability of each signal giv en the data as the Gaussian 
(iBond fc Critteudeull200l 

lnP(A x |A) = -i(A x - (A x |A)) t (C x )- 1 (A x - (A X |A)) 

-^TrlnC* -AHnV^, (6) 

where A are the observations, A x is the map vector of the 
component X, C x is the mode to mode correlation ma- 
trix for the component anc^TV is the number of modes in 
the map. The mean (A X |A) defines the Wiener-filtered 
map of the component given the observations 

A x = c x WA? ( 7 ) 

with variance about the mean 

(£A x £A xt |A) = C x - C X WC X . (8) 

In our case, the maps A are column vectors contain- 
ing the gridded visibility estimators (Paper IV), and 
sky plane maps can be obtained by Fast Fourier Trans- 
forming these estimators in the iw-plane. The corre- 
lation structure of the estimator grids are a necessary 
by-product of our gridding method and power spectrum 
estimation pipeline. The calculation of the correlations 
includes all the relevant information on the sampling 
structure and convolution of the fine-grained observa- 
tion plane. Eq. [HI therefore includes all aspects of the 
uncertainties in the Wiener-filtered map resulting from 
the observations. 
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Right Ascension (arcmin) 

Fig. 8. — The image produced by gridding a simulated observa- 
tion of a point source placed at the center of the 08 h field. 

Due to the complex weighting applied to the observed 
visibilities in our gridding scheme it is difficult to assign 
units to the maps produced by simply Fourier transform- 
ing the gridded ^v-plane. In order to obtain images with 
approximate normalizations in mJy/pixels, we reweight 
by dividing the vector at each lattice site by a second 
vector obtained by gridding a 1 mJy point source placed 
at the center of each field. It is important to note that 
this deconvolution is carried out in the coarse grained 
lattice and is not a full deconvolution of the effect of the 
dirty beam since no extra information is added to com- 
plete the ^i^plane when gridding. In Fig. [HI we show the 
image obtained by Fourier transforming the gridded uv- 
plane of a simulated observation of a point source placed 
at the center of the 08 h field. 

We show an example of the use of such filters in Fig. [31 
The sequence shows the result of applying different fil- 
ters to the 08 h deep field observation. The original image 
(upper left) is first filtered to obtain a map of the to- 
tal signal contributions which includes CMB and point 
sources (upper right). The CMB power spectrum ob- 
tained from the joint deep field analysis is then used as 
a template to obtain a map of the CMB contribution 
(lower left). The template includes the excess power 
above i ~ 2000. The post-subtraction residuals in the 
known point sources can also be separated into an im- 
age (lower right) which shows the total contribution of 
the modes which are projected out when estimating the 
power spectrum. The mean Wiener-filtered maps make 
no statement on the significance of any feature by them- 
selves and considering only the mean can be misleading 
since there is no information on the allowed fluctuations 
around the mean in the map itself. This is particularly so 
in the case of interferometric observations where the na- 
ture of the noise and sampling uncertainties means that, 
although the power due to single features is conserved 
their structure in the sky plane is complicated by the 
extended nature of the correlations in the uncertainties. 
As shown above, however, in the Gaussian limit we can 
assign confidence limits on any features using Eq. [51 

A useful method of visualizing the significance of the 
features is by creating constrained realizations of the flue- 
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Fig. 9. — Filtered images of the CBI 08 h deep field. The total signal plus noise image (top left) displayed as a 70' by 70 ; field. The 
white (dashed) circle shows the 45' FWHM of the primary beam. A total signal image (top right) is obtained by using the sum of CMB 
and point source foreground correlations as a template. The amplitude of the residual source background is set to a flux of 0.08Jy 2 sr _1 
in the filter. Using the joint deep field bandpowers as amplitudes in the filters we obtain the optimal image for the CMB (bottom left). 
This image encompasses the power attributed to the high-^ excess. The known point source residuals can also be filtered out of the data 
(bottom right). Colour scales are kept constant for the filtered images to ease comparison of the relative amplitudes of the components. 
The images are deconvolved by the response of a simulated source placed at the center of the field. 



tuations and comparing these to the mean. For exam- 
ple we can obtain a random realization of the fluctu- 
ations by taking Af = (C x - C X WC X ) 1 / 2 ^ where f 
are random Gaussian variates with unit norm. Adding 
these maps to the mean we get an idea of the allowed 
fluctuations about the mean. Strongly constrained fea- 
tures will be relatively unaffected by the fluctuations 
and represent the high signal-to-noise ratio limit where 
C x — C X WC X — > 0. Features measured at levels com- 
parable to the 'generalized' noise of the observations will 
be washed out by the fluctuation levels and represent the 
limit C x - C X WC X -> C x . The full pixel-pixel covari- 
ance matrix of the component map can be calculated by 
direct FFT of the covariance matrix or by taking Monte 



Carlo ensemble averages of the fluctuation realizations. 

The level of fluctuations allowed around the mean for 
the CMB component of Fig. El is shown in Fig. ^21 The 
four panel figure shows the Wiener filter of the compo- 
nent and the same map with three random realizations 
of the fluctuations added to it. As expected the CMB 
component is well constrained on the scales shown in the 
image since the signal-to-noise ratio is high: i.e. the 
main features in the map are stable with respect to the 
fluctuations. The features are therefore expected to be 
closely related to 'real' structures on the sky; but the 
nature of the observations, in particular the side-lobes 
of the point spread function (synthesized beam), means 
that the detailed structure of the features is not well 
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Fig. 10. — Another way to illustrate whether the features seen in the Wiener-filtered images are robust is to show a few fluctuations about 
the mean (Wiener) map, here for the 08 h deep field. The top left panel reproduces the mean CMB component map of Fig-EI The remaining 
three panels show the addition of three constrained realizations of the allowed fluctuations to the mean map. At the signal-to-noise ratios 
of the larger scale CMB observations the features are robust to the fluctuations. Note that these images are 'dirty maps': no attempt has 
been made to compensate for gaps in the m;- cover age. The features are a combination of main-lobe and side-lobe responses to the CMB. 



constrained and most likely does not reflect the precise 
shape of the objects. The filtering of the large scale power 
by the interferometric observations also makes the cor- 
respondence with the large scale features on the sky not 
as intuitive as with sky plane based measurements. 

Using the specific SZ models we introduced in §0]we 
can now extend the image analysis by template-filtering 
the images specifically for SZ contributions. As described 
above the only components modeled in the total correla- 
tion matrix are the CMB, instrumental noise and source 
components. We therefore construct the SZ template as 

/^pCMB 

cS z = Ec| z<^_. „, 

B 

where the amplitudes C^ z are obtained by filtering the 
average SZ power spectra by the deep field band window 
functions (see Paper II). 

The resulting SZ-filtered images should reproduce the 
structure in input SZ maps and we find that this is in- 
deed the case. The sequence of images in Fig. is an 
example of the filtering for one of the SPH maps (top 
row) and one of the MMH maps (bottom row). The left- 
hand panels show the input SZ maps (lead— trail) and 
the middle ones show the result of filtering mock data 
sets where the noise has been reduced to negligible lev- 
els. These provide a control set to compare with the re- 
constructed noisy observations shown on the right. The 
strongest features in both images reproduce the strongest 



SZ signatures in the input templates, although one has 
to take into account the tapering effect induced by the 
primary beam when comparing the input maps to the 
filtered data (see Paper IV). The tapering suppresses the 
amplitude of significant features further out in the beam 
with respect to those close to the center of the field. The 
noise level is such that the significance of the 'detections' 
is not obvious. These can be quantified by comparing 
with the allowed fluctuations about the mean Wiener fil- 
ter (Eq.®. 

We have also applied the SZ Wiener filter to the deep 
field data itself. These SZ-filtered images scale linearly 
with the overall normalization of the SZ power spectrum. 
Hence, by calibrating the SZ template filters with the ob- 
served deep high-^ power, we can obtain an image with 
the same signal-to-noise ratio as the observations. In 
fig El we show the effect of the SZ filter on the data split 
in two roughly equal halves for the 08 h and 20 h deep field 
observations. When comparing the two halves the noise 
will be uncorrelated, so any commmon stationary sig- 
nal of sufficient amplitude should appear in both halves. 
One sees no obvious strong signal-dominated stationary 
features, but there is one in each field, albeit with ampli- 
tudes comparable to the noise fluctuations. We conclude 
that the images are noise-dominated and there are no 
obvious features that could single-handedly account for 
the observed excess. This is in agreement with the ex- 
cess being a statistical measurement as opposed to a high 
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Fig. 11. — Simulated observations of SZ signals. A single SZ realization is shown for the SPH (top left) and MMH (bottom left) 
simulations. The area represents a (lead-trail) differenced field used as the signal in a simulated observation of the 08 h deep field. The 
maps are used to generate mock data sets with exact replications of the wu-coverage as the real data which are then filtered using the 
respective SZ templates. Mocks with close to negligible noise levels (top and bottom center) show the direct representation of the CBI 
processed input maps. The same realization but with the same noise levels as the observed data is in the top and bottom right panels. 
To within the accuracy allowed by the noise, the filtered images successfully reproduce the features of the input maps within the primary 
beam area. 



signal-to-noise observation. 

6. DISCUSSION 

We have presented estimates for the possible contri- 
bution from the SZ effect to the CBI deep field obser- 
vations which show evidence of an excess in power over 
standard primary CMB scenarios above £ ~ 2000. Our 
numerical simulations show that an amplitude of as ~ 1 
appears to give, on average, enough power in the SZ sig- 
nature to explain the excess. Current primary CMB and 
LSS data favor models with somewhat lower normaliza- 
tions: for flat, HST-ft and LSS priors, <j 8 = 0.88±g;g| 

for the June 2002 dataset and a 8 = 0.86±g;g| for the 
March 2003 compilation. Recent weak-lensi ng-on fy re- 
sults, e.g., 0.83 ±0.07 (|van Waerbeke et aljEg), are 
also quite compatible with CMB-only values. 

In assessing whether as ~ 0.9 is too low to explain the 
anomaly, the no n- Gaussian nature of the sample variance 
of the SZE should be included. T his issue has been par- 
tially addressed in the literature (|Coorav fc Shethll 2002: 
[ Springe!. White fc HerrinuistJ 120011 iKomatsu fc Seliakl 



2002). Goldstein etTal] 



have shown that a joint 



analysis of our CBI deep data with recent data from 
ACB AR (|Kuo et al.1 I2004D and BIMA {Dawson et all 
2002) which simultaneously fits the amplitudes of the pri- 
mary and SZ power spectra and roughly accounts for the 



non-Gaussianity of the SZE requires as = 0.96io'.?6? w ^ n 
a maximum likelihood at 1.00, consistent with the con- 
straints on erg from prima ry CMB and LSS data, as illus- 
trated in Figures [I] and Readhead et al. (2004) applied 
a similar joint analysis to that in Goldstei n et alJ (2003) 
to the bandpowers of the ACBAR, BIMA, and CBI 2- 
year data. An attempt to bracket the effects of non- 
Gaussianity was made by allowing the sample- variance 
component of the error bar to increase by a factor of up 
to four: the estimate of 0% hardly changed, going from 
0.98±g;g? to °- 98 ^o.08- IGoldstein et"all (I200I used a fac- 
tor of three to estimate the non- Gaussian effect. A bigger 
impact is found at the 2a level: with a uniform prior in 
as there is a tail skewed to low as, primarily because the 
amplitude being fit by the analysis is cr|, and the dis- 
tribution is roughly G aussian in that v a riable , as shown 
explicitly in Fig. 12 of iReadhead et al.1 (|2004l) . We note 
that the errors are predominantly from noise and not 
sample-variance for the CBI data, which is why varia- 
tion b y a factor of fou r does not have a huge effec t. 

The lGoldstein et al.1 ()2003ft and lReadhead et al.1 (2004) 
analyses took into account the competing effects of the 
damping tail of the primary CMB fluctuations and the 
rising SZ power spectrum over the entire £ range, and not 
just in the £ > 2000 regime where the "excess" is clearly 
seen. The assumption is that the power spectra are de- 




Fig. 12. — The split deep data Wiener-filtered with the SZ template, for the 08 h field (top) and the 20 h field (bottom). There is no clear 
evidence of a high signal-to-noise stationary feature. 



scribed by an offset-log-normal distribution, which has 
been shown to hold well for Gaussian s ignals and noise 
(|Bond. Jaffe fc KnoxEOQCl ISievers et aUfcoOfl . The dis- 
tribution becomes a Gaussian one in the noise-dominated 
regime and a log-normal one in the signal-dominated 
regime. If the signal comes from clusters and groups, the 
latter will not be accurate, and further work is needed 
to properly assess the limits derived from the "excess". 
What is needed is a full suite of simulations, guided by 
semi-analytic theory, to Monte Carlo the effects of non- 
Gaussianity rather than rely on expanded signal-errors. 

The excellent agreement among the runs using the two 
separate hydrodynamic algorithms is encouraging and 
lends support to the conclusions on the required nor- 
malization. We note that the simulations employed in 
this work were developed and run in dependently from 
each other and the CBI an alysis effort (B ond et al.l 2002: 
Zha ng. Pen fc Wangl 120021) . However, both simulations 



were adiabatic, with entropy increases occurring only by 
shock heating. The influence of cooling and entropy in- 
jection from supernovae, etc., on the predicted SZ spec- 
trum remains to be explored, as does the effect of such 
technical issues as numerical convergence as the simula- 
tion resolution is changed. 

The limited number of mock simulations we have car- 
ried out show that our experimental setup does repro- 
duce correctly the signal of simulated SZ foregrounds. 
This suggests the Gaussian assumption implicit in our 
bandpower estimation procedures does not break down 
at the SZ signal-to- noise ratios of the deep field observa- 
tions. 

If the excess is not due to the SZE, what could it be? 
Other possible origins of the excess were addressed in 
Paper II. It is inconsistent with adiabatic inflationary 
predictions for primary anisotropics on the damping tail 
only at the ~ 3a level. The excess clearly needs to be 
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further explored, both within the CBI dataset itself, and 
also by correlating with other observations of these deep 
fields. The statistical nature of the detected excess means 
we cannot identify single features responsible for the ex- 
cess in the deep fields even after Wiener-filtering the data 
with SZ derived templates. 

Our principal conclusion is that the CBI excess could 
be a result of the SZ effect for the class of ACDM con- 
cordance models if erg is in the upper range of values al- 
lowed by current CMB data. The simple Cg oc aKflbh) 2 
scaling (eq. |2j) shows that the lower <jg values preferred 
by the CMB data imply SZ signals below, but not too 
far below, the sensitivity obtained by the CBI deep ob- 
servations. This work also highlights how the signature 
of the SZ effect has great potential for constraining the 
amplitude <r 8 . The sensitive scaling of the CMB power 
with <jg will result in significant constraints even with 
large errors on the observed bandpowers. This further 



underlines the significance of 'blank field' observations at 
£ > 2000 which should reveal the SZ structure that neces- 
sarily lurks as a byproduct of CMB-normalized structure 
formation models. 
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